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^ ' Abstract 

This paper presents the first implementation of a two-way relay network based on 
the principle of physical-layer network coding. To date, only a simplified version of 
physical-layer network coding (PNC) method, called analog network coding (ANC), 
has been successfully implemented. The advantage of ANC is that it is simple to 
implement; the disadvantage, on the other hand, is that the relay amplifies the noise 
along with the signal before forwarding the signal. PNC systems in which the relay 
performs XOR or other denoising PNC mappings of the received signal have the 
^ I potential for significantly better performance. However, the implementation of such 

PNC systems poses many challenges. For example, the relay must be able to deal with 
symbol and carrier-phase asynchronies of the simultaneous signals received from the 
>• . two end nodes, and the relay must perform channel estimation before detecting the 

^ ' signals. We investigate a PNC implementation in the frequency domain, referred to as 



in 
o 



; FPNC, to tackle these challenges. FPNC is based on OFDM. In FPNC, XOR mapping 

^ ' is performed on the OFDM samples in each subcarrier rather than on the samples 

in the time domain. We implement FPNC on the universal soft radio peripheral 
(USRP) platform. Our implementation requires only moderate modifications of the 
packet preamble design of 802.1 la/g OFDM PHY. With the help of the cyclic prefix 
(CP) in OFDM, symbol asynchrony and the multi-path fading effects can be dealt 
^ I with in a similar fashion. Our experimental results show that symbol-synchronous 

^ ■ and symbol-asynchronous FPNC have essentially the same BER performance, for 

■ ■ ■ both channel-coded and unchannel-coded FPNC. 
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1. Introduction 

In this paper, we present the first implementation of physical-layer network coding 
(PNC) on the software radio platform. We believe this prototyping effort moves the 
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concept of PNC a step toward reality. Our implementation work also exposes and 
raises some interesting issues for further research. 

PNC, first proposed in jH, is a subfield of network coding j2| that is attracting 
much attention recently. The simplest system in which PNC can be applied is the two- 
way relay channel (TWRC), in which two end nodes A and B exchange information 
with the help of a relay node R in the middle, as illustrated in Fig. [TJ Compared 
with the conventional relay system, PNC could double the throughput of TWRC by 
reducing the needed time slots for the exchange of one packet from four to two ^ . In 
PNC, in the first time slot, end nodes A and B send signals simultaneously to relay R; 
in the second phase, relay R processes the superimposed signals and maps them to a 
network-coded packet for broadcast back to the end nodes. From the network-coded 
packet, each end node then makes use of its self information to extract the packet 
from the other end node [H, Is], 0]. 

Prior to this paper, only a simplified version of PNC, called analog network coding 
(ANC) jij , has been successfully implemented. The advantage of ANC is that it is 
simple to implement; the disadvantage, on the other hand, is that the relay amplifies 
the noise along with the signal before forwarding the signal, causing error propagation. 

To our best knowledge, the implementation of the original PNC based on XOR 
mapping as in [l| has not been demonstrated, even though it could have significantly 
better performance. A reason is that the implementation of XOR PNC poses a 
number of challenges. For example, the relay must be able to deal with symbol 
and carrier-phase asynchronies of the simultaneous signals received from the two end 
nodes, and the relay must perform channel estimation before detecting the signals. 

This paper presents a PNC implementation in the frequency domain, referred to as 
FPNC, to tackle these challenges. In particular, FPNC is based on OFDM, and XOR 
mapping is performed on OFDM samples in each subcarrier rather than the samples 
in the time domain. We implement FPNC on the universal soft radio peripheral 
(USRP) platform. Our implementation requires only moderate modifications of the 
packet preamble design of 802.11a/g OFDM PHY. With the help the cyclic prefix 
(CP) in OFDM, symbol asynchrony and the multi-path fading effects can be dealt 
with in a similar fashion. Our experimental results show that symbol-synchronous 
and symbol-asynchronous FPNC have nearly the same BER performance, for both 
channel-coded and unchannel-coded FPNC. 
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Figure 1: System model for physical-layer network coding. 
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Challenges 

In the following, we briefly overview the challenges of PNC, and the implementa- 
tion approaches taken by us to tackle them: 

Asynchrony 

There are two possible implementations for PNC: synchronous PNC and asyn- 
chronous PNC. In synchronous PNC, end nodes A and B have the uplink channel 
state information (CSI). They perform precoding and synchronize their transmissions 
so that their signals arrive at relay R with symbol and carrier-phase alignments. For 
high-speed transmission, such tight synchronization is challenging; in addition, timely 
collection of CSI is difficult in fast fading scenarios. 

Asynchronous PNC is less demanding. It does not require the two end nodes to 
tightly synchronize and precode their transmissions. In particular, knowledge of the 
uplink CSI is not needed at the two end nodes. The simplicity at the end nodes 
comes with a cost. Without precoding and synchronization of the two end nodes, 
their signals may arrive at the relay with symbol and carrier-phase misalignments. 
A key issue in asynchronous PNC is how to deal with such signal asynchrony at the 
relay 0,3. 

This paper focuses on the implementation of asynchronous PNC. To deal with 
asynchrony, our FPNC implementation makes use of OFDM to lengthen the symbol 
duration within each subcarrier. Then, independent XOR PNC mapping is performed 
within each subcarrier. OFDM splits a high-rate data stream into a number of lower- 
rate streams over a number of subcarriers. Thanks to the larger symbol duration 
within each subcarrier, the relative amount of dispersion caused by the multipath 
delay spread is decreased. The OFDM symbols of the two end nodes become more 
aligned with respect to the total symbol duration, as illustrated in Fig. [2J In par- 
ticular, if the relative symbol delay is within the length of the CP, the time-domain 
misaligned samples will become aligned in the frequency domain after DFT is applied. 
This property will be further elaborated later in Section 2. 

Channel Estimation 

For good performance of asynchronous PNC, the relay must have the knowledge 
of the uplink CSI. This has been the assumption in many prior works on PNC (e.g., 
[ll,[8|). This means that in implementation, the relay will need to estimate the channel 
gains. Most channel estimation techniques for the OFDM system assume point-to- 
point communication in which only one channel needs to be estimated. In PNC, the 
relay needs to estimate two channels based on simultaneous reception of signals (and 
preambles) from the two end nodes. This poses the following two problems in PNC 
that do not exist in point-to-point communication: 

• Channel estimation in a point-to-point OFDM system (e.g., 802.11 jol) is gener- 
ally facilitated by training symbols and pilots in the transmitted signal. If used 
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Figure 2: PNC with time asynchrony: (a) frequency-domain physical-layer network coding (FPNC); 
(b) time-domain physical-layer network coding (TPNC). 



unaltered in the PNC system, the training symbols and pilots from the two end 
nodes may overlap at the relay, complicating the task of channel estimation. 
In our implementation, we solve this problem by assigning orthogonal training 
symbols and pilots to the end nodes. The details will be given in Section 4. 

• It is well known that carrier frequency offset (CFO) between the transmitter 
and the receiver can cause inter-sub carrier interference (ICI) if left uncorrected. 
In a point-to-point system, CFO can be estimated and compensated for. In 
PNC, we have two CFOs at the relay, one with respect to each end node. 
Even if the two CFOs can be estimated perfectly, their effects cannot be both 
compensated for totally; the total elimination of the ICI of one end node will 
inevitably lead to a larger ICI for the other end node. To strike a balance, our 
solution is to compensate for the mean of the two CFOs (i.e., compensate for 
(CFO^ - CF0ij)/2). The details will be elaborated in Section 4. 

Joint Channel Decoding and Network Coding 

For reliable communication in a practical PNC system, channel coding needs to 
be incorporated. This paper considers link-by-link channel-coded PNC, in which the 
relay maps the overlapped channel-coded symbols of the two end nodes (1, to 
the XOR of the source symbols 0; after that, the relay channel-encodes the XOR 
source symbols to channel-coded symbols for forwarding to the end nodes. Such a 
link-by-link channel-coded PNC system has better performance than an end-to-end 
channel-coded PNC system 0, • 



"'^This process is called Channel-decoding-Network-Coding (CNC) in [10| because it does two 
things: channel decoding and network coding. Unlike the traditional multiuser detection (MUD) 
in which the goal is to recover the individual source information from the two end nodes, CNC 
aims to recover the XOR of the source information during the channel decoding process. CNC is a 
component in link- by-link channel-coded PNC critical for its performance 0, . 
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In our FPNC design, we adopt the convolution code as defined in the 802.11 
a/g standard. The relay first maps the overlapped channel-coded symbols to their 
XOR on a symbol-by- symbol basis. After that it cleans up the noise by (i) channel- 
decoding the XOR channel-coded symbols to the XOR source symbols, and then (ii) 
re- channel- CO ding the XOR source symbols to the XOR channel-coded symbols for 
forwarding to the two end nodes. 

The remainder of this paper is organized as follows: Section 2 details the delay 
asynchrony model of this paper. Section 3 presents the FPNC frame format design. 
Section 4 addresses the key implementation challenges. Experimental results are given 
in Section 5. Finally, Section 6 concludes this paper. 

2. Effect of Delay Asynchrony in Frequency Domain 

In asynchronous PNC, symbols of the two end nodes may arrive at the relay 
misaligned. We mentioned in the introduction that if the relative symbol delay is 
within the length of the CP in FPNC, then the time-domain misaligned samples 
will become aligned in the frequency domain after DFT is applied. This section is 
devoted to the mathematical derivation of this result. Here, we will derive a more 
general result that takes into account multi-path channels as well. 

2.1. Effective Discrete-time Channel Gains 

We consider the following multi-path channel model. Suppose that there are 
paths from node A to relay R with delays < < • • ■ < t^^'*"^ and corresponding 
channel gains a^, a^, . . . , The channel impulse response of A is gA{t) = 

Ma-I 

^ a\6{t — r^). Similarly, there are Mb paths from node B to relay R with delays 

1=0 

< < ■ ■ ■ < T^^^^^ and channel gains a]^, . . . , a^'^"^, with channel impulse 
response gB^t) = c^B^it — t^b)- Without loss of generality, we assume that frame 
A arrives earlier than frame B: specifically, 

rA<TB- Note that our model allows for 
the case where nodes A and B do not exactly transmit at the same time. If one node 
transmits slightly later than the other, we could simply add the lag time to all the 
path delays of that node. We assume that the net effect is such that the signal of 
A arrives earlier than the signal of -B, whether this is due to earlier transmission or 
shorter path delay of A . 

We first derive the effective discrete-time channel gains for the uplink in FPNC. As 
shown in Fig. [3l the discrete-time channel gains capture not just the continuous-time 
channel gains, but also the operations performed by pulse shaping and matched- 
filtering-and-sampling. Let us assume that the pulse shaping function p{t) is of finite 
length: specifically, we assume p{t) = for t < and t > Tp. The continuous-time 
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Figure 3: (a) Continuous-time channel model for PNC, in which and are the time 

domain source samples; 2/[n] is the time domain received samples; gA(t) and gsit) are the wireless 
multipath channel gains; p{t) is the pulse shaping function; w{t) is the receiver noise; and MF is the 
matched filter and sampler at the relay node, (b) Equivalent discrete-time channel model for PNC, 
in which /i^M and denote the equivalent discrete time channel impose response (i.e., effective 

discrete time channel gains), and w[n\ is the equivalent discrete-time noise term. 



baseband signal fed into the continuous-time channel is XA{t) = XA[n]p{t — riT). 

n=—oo 

The time domain received signal if only node A transmits is 

oo Ma-1 

VAit) = XA{t) * gA{t) + w{t) = Yl ^>A[n]p{t -t\- nT) + w{t), (1) 

n=-~oo 't=0 

where is the noise, assumed to be AWGN. Matched-Filtering (MF) and sampling 
are then performed on ([T]), by sampling at the first multipath channel tap of the uplink 
channel between node A and relay i?, to get the received samples 

oo 

yA[m] = j y{t)p{t -t\ + Tp- mT)dt 

— oo 

oo ( J Ma-1 1 

= ^ xaM < / Y'^AP{'^-A-nT)p{t-Tl + Tp-mT)dt\+w[m] 

n=-oo i=0 J 

oo 

= XA[n]hA[m — n] + w[m], (2) 



where w[m] = J w{t)p{t — r'^ + Tp — mT)dt. We see that the effective discrete-time 

— oo 

oo A/^ — 1 

channel of A is such that hAim—n] = f J2 (^APi't — Ta~ nT)p(t — t\ + Tp — mT)dt. 

— oo i=0 

Note that p{t-T\-nT)p{t-Tl+Tp-mT) = if \t\ -T\ + Tp~{m- n)T\ > Tp. 
In other words, hAlTn — n] =0 for (m — n)T > t^^~'^ ~ ''"a + and (m — n)T < 0. 
Define Da = [(r^^"' - r « + 2Tp) /T] . 

Let us now consider what if both end nodes transmit. The received signal at the 
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relay node is 



yR{t) = XA{t) * gA{t) + XB{t) * gsit) + w{t). 



(3) 



Sticking to the above MF that is defined with respect to first path delay oi A, we 
have 



y^[n] = XA[n] * hA[n] + Xsin] * hB[n] + w[n], 



(4) 



where /i^M = forn < and n > Da = \{rA'^~^ -^a + '^Tp) /T] , and hB[n] = 
for n < [(rO - + 2Tp) /T] and n>DB= [(r^ - + 2Tp) /T] . 



12 3 



Da-1 Da Da+1 



KM 



I 1 



12 3 Db-1 Db Db+1 

Figure 4: Example of delay spread in FPNC. 



2. 2. Delay- Spread- Within- CP Requirement 

The delay spread of node A is Da, and the delay spread of node B, with respect to 
time n = 0, is Db- We define the delay spread of the PNC system as (i.e., it combines 
the delay spreads of A and B into a potentially larger delay spread, as illustrated in 
Fig. DelaySpread) 



delay spread = max [Da, Db] ■ 



(5) 



The above derivation is general and do not have any requirement on the modulation. 
This subsection, we will present the OFDM modulated PNC system. In particular, 
we will present the "Delay-Spread- Within-CP Requirement" for FPNC. That is, we 
combine the Cyclic Prefix (CP) and Discrete Fourier Transform (DFT) to show that 
the time-domain symbol asynchrony of FPNC disappears in the frequency domain, 
when the uplink frames satisfy the Delay-Spread- Within-CP Requirement. 

First, let H[k] be the iV-point DFT of h[n] given by [lH (Note that, in the following 
derivation, we assume the subcarrier indices start from 0, i.e, k = 0, - ■ ■ ,N — 1) 



N-l 



H[k] = DFT{h[n]} = 



n\e 



0< k< N -1. 



(6) 



n=0 
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The A^-point circular convolution of Xn and hn is written as 

N-l 



y[n] = x[n] ®n h[n] = h[k]x[n — k]^, (7) 



k=0 



where [n — k]N denotes [n — k] modulo N. In other words, x[n — k]N is a periodic 
version oi x[n — k] with period A^. From the definition of DFT, circular convolution 



in time leads to multiplication in the frequency [11 



DFT{x[n]®Nh[n]} = X[k]H[k], < k < N - 1. (8) 

The channel output, as in (jlj), however, is not a circular convolution but a linear 
convolution. The linear convolution between the channel input and impulse response 
can be turned into a circular convolution by adding a special prefix to the input called 
a cyclic prefix (CP) fl2| . 

For FPNC, let H^lk] and Hslk] denote the frequency responses of the discrete- 
time channels, and let C denote the length of the CP. One OFDM symbol duration is 
then N + C. The CP for x^[n] is defined as xa[N — C], . . . , xa[N — 1]: it consists of 
the last C values of the xa [n] sequence. For each input sequence of length N, these last 
C samples are appended to the beginning of the sequence. This yields a new sequence 
x^^^^[n], -C<n<N-l,oi length N + C, where x^^^^[-C], . . . , a:2™^[iV - 
1] = xa[N - C], . . . ,xa[N - 1],xa[0], . . . ,xa[N - 1]. Note that with this definition, 
^OFD A/ ^ xaItiIn for —C < n < N — \, which implies that x^^^^^[n — k] = 
XA[n — k]N for —C <n — k<N — 1. 

Suppose Xa^'^^Iu] and x^^^^[n] are inputs to a discrete-time channel with im- 
pulse response hA[n] and hB[n], respectively. The channel output ynln], < n < N—1 
is then (assuming that the delay spread of FPNC max [Da, Db] is no large than the 
CP length C) 

yn[n] = x^^^^^n] * /i^[n] + xg^^^^n] * hsln] + w[n] 
c-1 c-1 
= hA[k]x'^^''^'[n -k] + Yl /iB[A;]a;g™[n - k] + w[n] 

k=0 k=0 
C-1 C-1 

= ^ hA[k]xA[n - k]N + ^ hB[k]xB[n - k]N + vj[n] 

k=0 k=0 

= XA[n\ ®Ar hA[n] + XB[n] (^n hB[n] + w[n], (9) 



where the third equality follows from the fact that for < A; < C — 1, x'^^^'^^[n — k] = 
XA[n — k]N for —C < n — k < N—1. Thus, by appending a CP to the channel 
input, the linear convolution associated with the channel impulse response ynln] for 
< n < — 1 becomes a circular convolution. Taking the DFT of the channel 
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output in the absense of noise then yields the following FPNC frequency domain 
digital expression: 

FPNC Frequency Domain Digital Expression: 

Define C as the length of the CP, and assuming FPNC delay spread = max [Da-, -Db] < 
C (where Da and Db are functions of the multipath delays t^^'^"^ and t^^~^, re- 
spectively) the received signal at subcarrier k is given by 

Y[k] = HA[k]XA[k] + HB[k]XB[k] + W[k], k = 0,...,N-l. (10) 

Note that the time-domain delay spread has been incorporated into HA[k] and HB[k] 
respectively. In FPNC, we will map Y[k] for each subcarrier k into the XOR, © 
Xb[/c]. This will be detailed in Section 4.3. The main point here is in (fTOj) . the signals 
of different subcarriers k are isolated from each other, and we only need to perform 
PNC mapping within each subcarrier. 

We remark that our discussion so far in this section has assumed the absence of 
CFO. When there is CFO, inter-carrier interference (ICI) may occur, and this will be 
further discussed in Section 4.1. 

3. FPNC Frame Format 
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Figure 5: FPNC preamble format. 

This section focuses on the PHY frame design to enable asynchronous operation, 
channel estimation, and frequency offset compensation in FPNC. As previously men- 
tioned, the asynchronous operation requires the PNC delay spread to be within CP. 
To ensure this, a simple MAC protocol as follows could be used to trigger near- 
simultaneous transmissions by the two end nodes. The relay could send a short 
polling frame (similar to the "beacon frame" in 802.11 that contains only 10 Bytes) 
to the end nodes. Upon receiving the polling frame, the end nodes then transmit. 
With this method, the symbols would arrive at the relay with a relative delay offset 
of \RTTa — RTTb\, where RTT is the round trip time, including the propagation 
delay and the processing time at the end nodes. This delay offset is not harmful to 
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our system as long as the sample misalignment of two end nodes is within the CP 
length. 

Given this loose synchronization, our training symbols and pilot designs described 
below can then be used to facilitate channel estimation and frequency offset compen- 
sation in FPNC. We modify the PHY preamble design of 802.1 la/g for FPNC. The 
overall FPNC frame format is shown in Fig. |5l The functions of the different com- 
ponents in the PHY preamble are described in the next few subsections. 

3.1. FPNC Short Training Symbol 

In 802.11, the short training symbol (STS) sequence contains 160 time-domain 
samples, in which 16 samples form one STS unit (sts) for a total of 10 identical 
units, as shown in Fig. El FPNC adopts the same STS sequence as in 802.11, as 
illustrated in Fig. |6l The STS sequence is used by the relay node to perform the 
sample timing recovery on the received frame. In particular, the relay node applies 
a cross-correlation to locate the sample boundary for the long training symbols that 
follow the STS sequence. The normalized cross- correlation is defined as follows: 



Z\n] 



L-l 

{sts*\i]yR[n + i]) 

i=0 



L-l 

E iyR[n + i]y*R['n + i]) 

i=0 



fill 



where n is the received sample index, yR[n] is the n-th sample at the relay R, and 
L = 16 is the length of each sts. For FPNC, this cross-correlation will result in 20 
peaks over the STS sequences (see Fig. [8]) of the two frames if the frames are not 
synchronized. From this profile of peaks, we can identify the last two peaks. If the 
Delay-Spread- Within-CP requirement is satisfied, then the last two peaks must be 
the last peaks of A and B, respectively. This is because the CP as well as the sts are 
of 16 samples in length. From there, we could locate the boundaries of the LTS of 
A and B that follow. Note that when the STS sequences of nodes A and B overlap 
exactly, we will have ten peaks only. In this case, the LTS boundaries of A and B also 
overlap exactly, and we simply use the last peak to identify the common boundary. 
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Figure 8: Cross-correlation of the STS for the uphnk of FPNC. 



3.2. FPNC Long Training Symbol 

With reference to Fig. [71 the 802.11 LTS sequence contains 160 time-domain 
samples in which there is a CP followed by two identical LTS units, Its. The receiver 
uses the LTS sequence to perform channel estimation and CFO compensation. 

For FPNC, in order to estimate two uplink channel gains, we design the LTS so 
that it contains twice the length of LTS in 802.11a/g, as shown in Fig. [3 In Fig. [Tj 
we intentionally show the case in which the LTS sequences of the two end nodes are 
not exactly synchronized. Note that we change the 802.11 LTS design by shortening 
its original CP length from 32 to 16 to make sure that the two Its units of B will not 
overlap with the data of A that follows under the condition that the delay spread is 
less than the CP length of 16. This does not impose additional requirement on the 
delay spread, since the CP of the data OFDM symbols in 802.11a/g (and FPNC) have 
only 16 samples anyway (i.e., the delay spread must be within 16 samples anyway). 
Section 4 will detail the CFO compensation and channel estimation methods for our 
implement at ion . 
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Figure 9: Long Training Symbol design for FPNC (time domain). 
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3.3. FPNC Pilot 

There are four pilots for each OFDM symbol in 802.11, as shown in Fig. [91 The 
four pilots are used to fine-tune the channel gains estimated from LTS. In a frame, 
there are multiple OFDM symbols, but only one LTS in the beginning. In practice, 
the channel condition may have changed by the time the later OFDM symbols arrive 
at the receiver. That is, the original channel gains as estimated by LTS may not be 
accurate anymore for the later OFDM symbols. The pilots are used to track such 
channel changes. 

In FPNC, we design the FPNC pilots of nodes A and B by nulling certain pilots 
to introduce orthogonality between them, as shown in Fig. O As will be detailed 
in Section 4.2, this allows us to track the channel gains of A and B separately in 
a disjoint manner in FPNC. We conducted some experiments for a point-to-point 
communication system using the two-pilot design rather than the four-pilot design. 
We find that for our linear interpolation channel tracking scheme described in Section 
4.2, the BER performances of the two-pilot and four-pilot designs are comparable for 
BPSK- and QPSK-modulated systems. 

4. Addressing Key Implementation Challenges in FPNC 

We next present our methods for carrier frequency offset compensation, channel 
estimation, and FPNC mapping, assuming the use of the PHY frame format presented 
in Section 3. 

4.1. FPNC Carrier Frequency Offset (CFG) Compensation 

For CFO compensation, we first estimate the two independent CFOs (namely 
CFOa and CFOb) caused by the carrier frequency offsets between nodes A and B 
and relay i?, respectively. We then compensate for the mean of the two CFOs (i.e., 
CFOpNc = {CFO A + CFOb) /2 ). The details are presented below. 

4.1.1. CFO Estimation 

For the uplink phase, when there are CFOs, the received frames at relay R will 
suffer from time-varying phase asynchronies. We need to compensate for the CFOs 
to alleviate inter-carrier interference (ICI) among data on different subcarriers. 

Recall that in Section 3, we mentioned that a loose synchronization MAC protocol 
can be used to ensure that the difference of the arrival times of the frames from nodes 
A and B are within CP. That means that the LTSs from nodes A and B will overlap 
with each other substantially, with the non-overlapping part smaller than CP (see 
Fig. [5]). Recall also that we introduce orthogonality between the LTSs of nodes A 
and B so that when the LTS units in A are active, the LTS in B are zeros, and vice 
versa, as shown in Fig. [71 This allows us to separately estimate CFO a and CFOb- 
Without loss of generality, in the following we focus the estimation of CFOa using 
LTSa. 
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CFOa is given by A/a = Ja — fn (i-e., the difference in the frequencies of the 
oscillators of node A and relay R). We define the normalized CFOa to be (pA = 
27rA/A^ , where T is the duration of one OFDM symbol, and N is the number of 
samples in one OFDM symbol not including CP. In other words, (pA is the additional 
phase advance introduced by the CFO from one sample to the next sample. 

To estimate (pA, we multiply one sample in the first unit of LTSa (see Fig. [7]) by 
the corresponding sample in the second unit of LTSa to obtain {y^^'^[Ti]) y^^^['n + 

N]. Then, angle (^{un^^'^in]) Vfi^^'^in + N]j E (—11,71) is given by 

angle [{yR^'^in])* y^/'^[n + iV]) + 2m7r = N(f)A, (12) 

where m e {. . . , -2, -1, 0, 1, 2, . . .}. 

For our experimental test-bed, USRP, we found that the accuracies of the onboard 
oscillators are such that they do not induce large CFOs so that m = (interested 



readers are referred to [ij] for CFO estimation when m 7^ 0.). Hence, we could write 
f lT^ as follows: 

NcPa = angle ((yf """M)* ^-^[n + iV]) . (13) 

Strictly speaking, (fT3!) is an expression for the noiseless case. Because of noise, 
angle (^{y^'^^^in]) y^^^^in + N]^ for different n G {0, . . . , — 1} could be different. 

Thus, in our computation, we first obtained 0a N = angle (y{y^^'^ [n\) y^^^ [n + N] 
for 77, = 0, . . . , — 1, and then estimate (pA by 



(pA = median [(t)A\'n\ ) . (14) 

ne{0,...,Af-l} V / 

We obtain (ps similarly. 

The reason we use the median CFO instead of the mean values is that we find 
the median is more stable. In particular, some samples of 0a i?^] are outliers that 
appear to be caused by unknown errors of significant magnitudes. We will show the 
BER results comparing the use of mean and median for CFO compensation (in Fig. 

mh))- 

4- 1-2. Compensation for Two CFOs 

In FPNC, we adopt the mean of the two CFOs for compensation purposes: 

a + ^b)/2. (15) 

Experimental results show (see Fig. [TOT a)) that compensation by the mean (p in 
f lT5]) is better than compensation by either 0a or (ps- We believe a theoretical study 
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(b) 



Figure 10: Comparison of tlic different CFO compensation methods: (a) CFOa v.s. 
(b) mean CFO v.s. median CFO. 



CFOa+CFOb . 

2 ■ 



to explore and compare different compensatation methods may be worthwhile in the 
future. As far as we know, there have been no theoretical treatments of compensating 
for two CFOs. Fig. [TOT b) shows the BER performance of using median for the 
estimate of cf) a oi (f)B as in ( IT4l) . versus using mean. It shows that he use of median 
results in better performance. 

After compensation, our received data in the time domain is given by 

VbH = yR[n]e-^''t (16) 

In the frequency domain, we have 

YR[k] = DFT{yn[n]). (17) 

We should emphasize that the computation complexity of FPNC CFO compensation 
is exactly the same as that of point-to-point communication, thus real-time decoding 
is possible. 

4-2. FPNC Channel Estimation 

In this subsection, we present the channel estimation and tracking method for 
FPNC. Note that CFO compensation was performed on the time-domain signal. For 
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channel estimation, however, we are interested in the channel gains for different sub- 
carriers in the frequency domain. This means that channel estimation will be per- 
formed after DFT. Thus, in the following we look at the signal after CP removal and 



For FPNC channel estimation, we use the LTS to obtain a first estimate. Pilots 
are used to obtain additional estimates for later OFDM symbols within the same 
frame. In the following, we consider channel estimation of HA[k]. Estimation of 
Hslk] is performed similarly. 

For channel estimation based on LTS, define one FPNC LTS unit of node A (i.e., 
with respect to Fig. [71 one unit is Its a) in the frequency domain as X^^[k], where 
k = 0, . . . , N — 1. Based on the first unit of Its a the received frequency domain 



LTSa (i.e., Y^'^^^ik] = DFT {y^^^^[n])), we perform channel estimation of HA[k] 



As mentioned in Section 3, each LTS contains two identical units in our design. 
The uplink channel gain HA[k] between node A and relay R is estimated by taking 
the average of the two units results 



In general, the channel may have changed from the first OFDM symbol to the last 
OFDM symbol within the same frame. The estimate based on LTS in f|T9|) applies 
only for the earlier symbols. Pilots are used to track the channel changes for later 
symbols. Our pilot design was shown in Fig. Pilot. In each FPNC OFDM symbol, 
there are two pilots per end node. Note from Fig. [Hlthat the two pilots of node A and 
the two pilots are node B are positioned at different subcarriers and non-overlapping 
in the frequency domain. Therefore, we could separately track the changes in HA[k] 
and Hslk]. In the following, we consider the tracking of HA[k]. Tracking of Hslk] 
can be done similarly. 

Let k' and k" denote the subcarriers occupied by the two pilots of A. Consider 
OFDM symbol m. Let Y^[k'] and Y^[k"] be the received signal in the frequency 
domain. Because the pilots of A and B do not overlap, Y^[k'] and Y^[k"] contain 
only signals related to the pilots of A. We first multiply YJ^[k'] and Y^[k"] by 
(i7J[A;'])-i and {H'^[k"])-^ obtained from ([19]), respectively. Let PaW] and PaW] 
be the two pilots. Then, we compute 



DFT. 



as follows: 



HA[k] 




(18) 




(19) 



AH^[k'] = m[k'])-'Ynk'] / PaWI 
AH'^lk"] = {H'^[k"])~'^Y^[k"] / PA[k"]. 



(20) 
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After that, we perform linear fitting to obtain AiJ™ [A;] for k ^ k', k" , as follows: 

' AH'X[k"] - AH^'[k'] 



AH'^ik] = AH'^ik'] + 



k" - k' 



{k-k'). 



(21) 



To obtain the final channel estimation for the m-th OFDM symbol, we compute 

H^[k] = H^^[k] ■ AR'^^ik]. (22) 
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Figure 11: Link-by-link channel-coded PNC, including channel-decoding and network coding (CNC) 
process and channel encoding. 



4.3. FPNC Mapping 

For reliable communication, channel coding should be used. Channel codin g in 
PNC systems can be either done on an end-to-end basis or a link-by-link basis |lOl . 
0|. The latter generally has better performance because the relay performs channel 
decoding to remove noise before forwarding the network-coded signal. 

The basic idea in link-by-link channel-coded PNC is shown in Fig. [11] It con- 
sists of two parts. Let denote the vector representing the overall channel-coded 
overlapped frames received by relay R. The operation performed by the first part is 
referred to as the Channel-decoding and Network- Co ding (CNC) process in 10|]. It 
maps Yr to Sa® Sb, where Sa and Sb are the vectors of source symbols from nodes 
A and B, respectively, and the © operation represents symbol- by-symbol XOR oper- 
ation across corresponding symbols in Sa and Sb- Note that the number of symbols 
in Yji is more than the number of symbols in Sa © Sb because of channel coding. 
Importantly, CNC involves both channel decoding and network coding. In particular, 
CNC channel- decodes the received signal Yr not to Sa and Sb individually, but to 
their XOR. The second part can be just any conventional channel coding operation 
that channel code Sa © Sb to Xji = C{Sa © Sb) for broadcast to nodes A and B, 
where C{*) is the channel coding operation. 



As mentioned in 10 



and 1^, the CNC component is unique to the PNC system, 
and different designs can have different performance and different implementation 
complexity. We refer the interested readers to jil for a discussion on different CNC 
designs. 

In this paper, we choose a design that is amenable to simple implementation, as 
shown in Fig. [T2j We refer to this CNC design as XOR-CD. In this design, any linear 
channel code can be used. In our implementation, we choose to use the convolutional 
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Figure 12: XOR-CD design for CNC. 



code. In XOR-CD, the channel-decoding and network coding operations in CNC are 
performed in a disjoint manner. As shown in Fig. [121 based on the CFO-compensated 



fc=0,l,. 



We then 



YfjJ/c] obtained as in (ITT]) , we obtain the overall vector Y/j = (^^[fc]) 
perform symbol-wise PNC mapping to get an estimate for the the channel-coded 



XOR vector Xa®Xb = (X^l/e] © ^^[/e]) 



fc=0,l, 



where Xa = {XA[k]) 



k=0,l,. 



and 



Xb = {XB[k]), 

=0 1 ... channel-coded vectors from A and B, respectively. We 

assume the same linear channel code is used at nodes A, B, and R. Note that since 
we adopt the convolutional code, C{*) is linear. Therefore, we have Xa © Xb = 
C{Sa) © C{Sb) = C{Sa © Sb), and thus the same Viterbi channel decoder as used 
in a conventional point-to-point communication link can be used in the second block 
of Fig. M 

The mapping in the first block in Fig. [12] could be performed as follows. Based 
on the channel gains estimated in ( 122|) . we could perform the XOR mapping for the 
k-th subcarrier in the m-th OFDM symbol (assuming BPSK modulation) according 
to the decision rule below: 



exp 



\Y]^[k]-H^[k]-H^[k]\ 



2a^ 



+ exp 



\Yg'[k]+H^[k]+H^[k]\ 



2a^ 



X^[k]=l 



exp 



\Y^[k]+Hl-[k]^H--[k]\ 



+ exp 



\Yl^[k]~H^[k]+WS[k]\ 



(23) 



where we have assumed Gaussian noise with variance o^. The computation complex- 
ity in (l23lll^. however, is large. In our implementation, we adopt a simple "log-max 
approximation" [l5[ (i.e., log(^. exp(2;j)) maxzj) that yields the following decision 

i 

rule: 



min - H^[k] - H^[kf , \Yl^[k] + H^[k] + H^[k]\^] 

^ min \ F™[A;] + H'X[k] - H%[k]\ , |Y^[A;] - H'Xlk] + H^^W 

X^[k] = l ^ 



(24) 



^Note that (US]) is similar to (7) in Ref. 

\H^[k]\^\RS[k]\ 



01, except that here we allow for the possibility that 
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This decision rule can also be interpreted as in Table [Tj where 

U= arg min {\Y;^[k] - Uf} . (25) 

Ue{±H^[k]±H]!}[k]} 



Tabic 1: XOR mapping with BPSK modulation in FPNC. 



U = arg min {|y^"[A;] - U\^} 

U&{±H^[k]±H^[k]} 




X^[k] = X^[k] (B X'Slk] 






1 






-1 


-Hm+H^[k] 




-1 


-H:^[k] - H^[k] 




1 



Note here that this decision rule could be used even for non-Gaussian noise. This 
is because ([25]) corresponds to finding the nearest point in the constellation map 
(constructed by combining the two end nodes' channel gains). 

Based on the XORed samples detected using the decision rule of Table [H we then 
perform the channel decoding to get the XORed source samples. In our implemen- 
tation, we use a Viterbi decoder with hard input and hard output. In general, a soft 



Viterbi algorithm could also be used for potentially better BER performance [16 



5. Experimental Results 

This section presents details of our FPNC implementation over the software radio 
platform and the experimental results. 

5.1. FPNC Implementation over Software Radio Platform 

We implement FPNC in a 3-node GNU Radio testbed, with Software Defined 
Radio (SDR). The topology is shown in Fig. [TJ Each node is a commodity PC 



connected to a USRP GNU radio 17 



Hardware: We use the Universal Software Radio Peripheral (USRP) [18| as our 
radio hardware. Specifically, we use the XCVR2450 daughterboard operating in 
the 2.4/5GHz range as our RF frontend. We use the USRPl motherboard for 
baseband data processing. The largest bandwidth that USRPl could support 
is 8MHz. In our experiment, we use only use half of the total bandwidth for 
FPNC (i.e., 4MHz bandwidth). 

Software: The software for baseband signal processing is based on the open 



source of GNURadio project [17|. We build our system by modifying the 802. llg 



transmitter implementation in the FTW project [19|. The FTW project [20 



however, does not have a 802. llg receiver. Therefore, we develop our own 
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OFDM receiver, designed specifically to tackle various issues in the FPNC sys- 
tem, such as CFO estimation and compensation, channel estimation, and CNC 
processing as presented in Section 4. 

5.2. Experimental Results 

We conduct our experiments over the channel one of 802.11g, with 2.412GHz being 
the central frequency. For each transmitter power level (we vary the SNR from 5dB 
to 20dB), we transmit 1000 packets and examine the resulting BER performance. 
Both the symbol-synchronous and symbol-asynchronous cases are investigated. The 
packet length is ISOOBytes, which is a normal Ethernet frame size. 

5.2.1. Time-Synchronous FPNC versus Time- Asynchronous FPNC 

In Section 2, we derived theoretically that as long as the Delay-Spread- Within-CP 
requirement is satisfied, FPNC will not have asynchrony in the frequency domain. Of 
interest is whether this reduces the asynchrony penalty in practice. In our first set 
of experiments, we investigate this issue. We study both unchannel-coded as well as 
channel-coded FPNC systems. 

To create different levels of time asynchrony, we adjust the positions of the end 
nodes. One of the set-ups corresponds to the perfectly synchronized case (the STS 
correlation has only ten peaks in the perfectly synchronized case: see Section 3). Fig. 
[TST a) shows the BER-SNR curves for the synchronous case, and Fig. [TST b) shows 
the curves for the asynchronous case with eight samples offset between the early and 
late frames. Note that this asynchrony still satisfies the Delay-Spread- Within-CP 
requirement because the CP has of 16 samples. We find that the performance results 
of the asynchronous cases with other time offset to be similar, and we therefore present 
the results of the eight-sample offset only. 

From Fig. [TST a) and (b), we see that the asynchronous FPNC has essentially the 
same BER performance as that of the synchronous FPNC. Hence, we conclude that 
FPNC is robust against time asynchrony as far as BER performance is concerned. 

5.2.2. FPNC versus Other Approaches for TWRC 

Our next set of experiments is geared toward the comparison of FPNC with other 
TWRC schemes. Recall that FPNC TWRC is a two-phase scheme using two time 
slots for the exchange of a pair of packets between two end nodes. We consider the 
following two additional approaches j3| 

• SNC: The straightforward network coding (SNC) scheme makes use of conven- 
tional network coding at the higher layer using three time slots. In SNC, node 
A transmits to relay R in the first time slot; node B transmits to relay R in 
the second time slot; relay R then XOR the two packets from A and B and 
transmits the XOR packets to nodes A and B in the third time slot. 
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Figure 13: BER of FPNC with and without sample synchronization. The 95% confidence intervals 
are marked in the figures. Note that the BER here is related to whether the XOR bit is decoded 
correctly, not whether the individual bits from the two end nodes are decoded correctly. 



• TS: Traditional scheduling (TS) scheme uses four time slots. In the first time 
slot, node A transmits to relay R; in the second time slot, relay R forwards the 
packet from A to node B. Similarly, the packet from node B to node A uses 
two additional times slots for its delivery. 

Our overall goal is the compare the throughputs of the three schemes. To derive 
the throughputs, we first measure the following three frame-error rates: 

1) FERpNc = p;^'^"'^'^^^: frame -error rate of the uplink of FPNC. 



2) FER, 



puphnk,SNC_^ frame-error rate computed from the two uplink time 



slots in SNC 

3) FERp2p = Pf^^'- frame-error rate of a point-to-point link. 

Channel coded systems are considered in our implementation. All three systems use 
the convolutional channel code with 1/2 coding rate, as specified in the 802.11a/g 
standard ^1. 
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Note for 1) and 2), FERpnc and FERsnc refers to the error rate of the XOR 
of the two source frames. That is the error rate for the frame Sr = Sa (B Sb- 
For FERsnc, we gather the decoded Sa and Sb from the two uphnk time slots, 
and the compute their XOR before checking whether there is an error in the XOR 
frame. For FERpnc, the CNC scheme as described in Section 4.3 is used to decode 
Sr directly based on the simultaneously received signals. In Fig. ITW a). we plot 
FERpNC: FERsnc^ and FERp2p versus SNR obtained from our experiments. The 
FER measurements are all from channel-coded systems. 

The throughputs per direction of the three TWRC schemes are computed as 
follows: 

rph ^ ^^ puplmk,PNC\f'i tdP2P\ 

^ IIPPNC = -^[i- - )[^ ~ -^f J' 

ThsNc = 1(1 - p;^'-'='^^^)(i - ppp), 

ThTs = l{l-Ppy- (26) 

In Fig. fW b). we plot the throughputs {Thppjqc^ ThcNC, and THts) of FPNC, 
SNC and TS versus SNR based on the FERppNc, FERcnc, and FERts in Fig. 
[TW a). With reference to Fig. [111(b), for the high SNR regime (above 19dB), the 
throughput of PNC is approximately 99% higher than that of the TS scheme, and 
49% higher than that of the SNC scheme. This is essentially the same as the ideal 
100% and 50% throughput gains derived by slot counting in [l| (i.e., the error- free 
case), with the difference that we have channel coding here to ensure reliable com- 
munication. If we use the guideline that the common decodable 802.11 link usually 
works at an SNR regime that is higher than 20dB 21, [sj, we can conclude that our 
FPNC implementation has very good performance in this regime. We note that for 
this regime, (sf mentions that ANC can achieve 70% and 30% throughput gains rel- 
ative to TS and SNC. Hence, FPNC has better performance in this SNR regime by 
comparison. 

We note from Fig. [TWb) that the performance of FPNC is not as good as that 
of SNC or TS at the lower SNR regime (say below 17dB). This is most likely due 
to our specific implementation of FPNC in this paper rather than a fundamental 
limitation of FPNC in general. In particular, recall that we implement the CNC 
function in FPNC mapping (see Section 4.3) using the so-called XOR-CD approach. 
In XOR-CD, (i) we first perform XOR mapping for the channel-coded symbol pairs 
from the two end nodes; (ii) after that channel decoding is applied on the channel- 
coded XOR symbols to get the XOR of the source symbols. Step (i) loses information 
that could be useful for the decoding of the XOR of the source symbols, and may 
lead to inferior performance in the low SNR regime. This phenomenon is explained in 
, and an joint CNC scheme [l3[ 0] for the PNC system can potentially achieve 



better performance than the XOR-CD scheme implemented in this paper, at the cost 
of higher implementation complexity. 
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Figure 14: Frame error rate and throughput comparison of FPNC with straight-forward network 
coding and traditional secluding, (a) FER comparison of three approaches; (b) throughput compar- 
ison of three approaches. 



6. Conclusion and Future Work 

This paper presents tlie first implementation of a PNC system in which the relay 
performs the XOR mapping on the simultaneously received signals as originally en- 
visioned in [l| . In particular, in our implementation, the XOR mapping is performed 
in the frequency domain of an OFDM PNC system. We refer to the OFDM PNC 
system as FPNC. The implementation of FPNC requires us to tackle a number of 
implementation challenges, including carrier frequency offset (CFO) compensation, 
channel estimation, and FPNC mapping. 

A major advantage of FPNC compared with PNC in the time domain is that 
FPNC can deal with the different arrival times of the signals from the two end nodes 
in a natural way. We show by theoretical derivation that if the simultaneously received 
signals in FPNC have a maximum delay spread that is less than the length of the 
OFDM cyclic prefix (CP), then after the Discrete Fourier Transform, the frequency- 
domain signals on the different subcarrier are isolated from each other. That is, in the 
frequency domain, the signals are synchronous. Then, straightforward XOR mapping 
can be applied on the different subcarrier signals separately in a disjoint manner. 
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To validate the advantage of FPNC, we present experimental results showing that 
time-domain symbol asynchrony does not cause performance degradation in FPNC. 

To date, most work related to PNC focuses on its potential superior performance 
as derived from theory. In this paper, we evaluate the throughput gain of PNC 
relative to other two-way relay schemes. Our implementation indicates that PNC 
can have a throughput gain of 99% compared with traditional scheduling (TS), and 
a 49% throughput gain compared with strait- forward network coding (SNC), in the 
high SNR regime (above 19 dB) in which practical technology such as Wi-Fi typically 
operates. 

Going forward, there are many rooms for improvement in our FPNC implemen- 
tation. In this paper, when faced with alternative design choices, we opt for imple- 
mentation simplicity than performance superiority. For example, we choose to use 
a simple PNC mapping method called XOR-CD in this paper, which is simple to 
implement but has inferior performance compared with other known methods j3] in 
the low SNR regime. In addition, our implementation exercise reveals a number of 
problems with no good theoretical solutions yet, and further theoretical analysis is 
needed; in such cases, we use simple heuristics to tackle the problems. For example, 
CFO compensation for FPNC is an area that is not well understood yet, because we 
have to deal with CFOs of more than one transmitter relative to the receiver. In this 
paper, we simply compensate for the mean of the CFOs of the two end nodes. Better 
methods await further theoretical studies. Last but not least, we base our design on 
the 802.11 standard to a large extent with only moderate modifications. If we do not 
limit our design within the framework of 802.11, there could be other alternatives 
with potentially better performance. 

References 

[1] S. Zhang, S. C. Liew, and P. P. Lam. "Hot Topic: Physical Layer Network 
Coding," in Proc. 12th MoMCom, pp. 358-365, NY, USA, 2006. 

[2] R. Ahlswede, N. Cai, S.-Y. R. Li, and R. W. Yeung, "Network Information 
Flow," IEEE Trans, on Information Theory, vol. 46, pp. 1204-1216, 2000. 

[3] S. Zhang, S. C. Liew, and L. Lu, "Physical Layer Network Coding Schemes over 
Finite and Infinite Fields," in Proc. IEEE GLOBECOM 2008. 

[4] S. C. Liew, S. Zhang, and L. Lu "Physical-Layer Network Coding: Tutorial, 
Survey, and Beyond," submitted to Phycom. 

[5] S. Katti, S. GoUakota, and D. Katabi, "Embracing Wireless Interference: Analog 
Network Coding," ACM SIGCOMM, 2007. 

[6] S. Zhang, S. C. Liew, and P. P. Lam, "On the Synchronization of Physical-Layer 
Network Coding," Information Theory Workshop, Oct 2006. 



23 



[7] L. Lu, S. C. Liew, and S. Zhang, "Optimal Decoding Algorithm for Asynchronous 
Physical-Layer Network Coding," in Proc. IEEE ICC 2011. 

Rossetto, F. and Zorzi, M., "On the Design of Practical Asynchronous Physical 
Layer Network Coding," IEEE 10th Workshop on SPAWC, 2009 

[9] IEEE 802.11-2007. Wireless LAN Medium Access Control (MAC) and Physical 
Layer (PHY) Specifications. 

[http : //standards . ieee . org/getieee802/download/802 . 11-2007 . pdf 



[10] S. Zhang, and S. C. Liew, "Channel Coding and Decoding in a Relay System 
Operated with Physical-Layer Network Coding," IEEE J on Selected Areas in 
Commun. vol. 27, no. 5, pp. 788-796, June 2009. 

[11] A. V. Oppenheim, et al., Discrete-Time Signal Processing (3rd Edition), Prentice 
Hall, 2009. 

[12] A. Goldsmith, et al.. Wireless Communications, Cambridge University Press, 
2005. 

[13] L. Lu, and S. C. Liew, "Asynchronous Physical-Layer Network Coding," 
|http : //arxiv . org/abs/1105 . 3144[ 

[14] S. G. Glisic, Advanced Wireless Networks: 4G Technologies (second eidtion), 
John Wiley & Sons, 2007. 

[15] S. Boyd, and L. Vandenberghe, Convex Optimization, Cambridge University 
Press, 2004. 

[16] B. Sklar, Digital Communications: Fundamentals and Applications (2nd Edi- 
tion), Prentice Hall PTR, 2003. 

[17] G. FSF. GNU Radio - Gnu FSF Project, 



http : / / gnuradio . org/ redmine/ wiki/gnuradio 



[18] Ettus Inc. Universal Software Radio Peripheral, 
|http : / /www ■ ettus . com/ 1 

[19] P. Fuxjaeger, et al., FTW IEEE802.11a/g/p OFDM Frame Encoder, 



|https : //www. cgran. org/wiki/f tw80211of dmtx 



[20] P. Fuxjer et al., "IEEE 802. lip Transmission Using GNURadio," in Proc. IEEE 
WSR, 2010. 

[21] J. Geier. SNR Cutoff Recommendations, 2005. 



http : //www. wi-f iplanet . com/tutorials/article .php/3468771 



24 



